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Abstract 

An extensive theoretical search for the proton magic number in the superheavy valley beyond Z =82 and correspond- 
ing neutron magic number after =126 is carried out. For this we scanned a wide range of elements Z - 1 12 - 130 
and their isotopes. The well established non-relativistic Skryme-Hartree-Fock and Relativistic Mean Field formalisms 
with various force parameters are used. Based on the calculated systematics of pairing gap, two neutron separation 
energy and the shell correction energy for these nuclei, we find Z = 1 20 as the next proton magic and N= 1 72, 182/1 84, 
208 and 258 the subsequent neutron magic numbers. 
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After the discovery of artificial transmutation of ele- 
ments by Sir Ernest Rutherford in 1919 01], the search 
for new elements is an important issue in nuclear sci- 
ence. The existence of elements beyond the last heavi- 
est naturally occurring ^■'^U, i.e., the discovery of Nep- 
tunium, Plutonium and other 14 elements (transuranium 
elements), which make a separate block in Mendeleev's 
periodic table was a revolution in the Nuclear Chem- 
istry. This enhancement in the periodic table raises a 
few questions in our mind: 

• Whether there is a limited number of elements that 
can co-exist either in nature or can be produced 
from artificial synthesis by using modern technique 

? 

• What is the maximum number of protons and neu- 
trons that of a nucleus ? 

• What is the next double shell closure nucleus be- 
yond 2"8pb ? 

To answer these questions, first we have to under- 
stand the agent which is responsible to rescue the nu- 
cleus against Coulomb repulsion. The obvious reply is 
the shell energy, which stabilises the nucleus against 
Coulomb disintegration f^. Many theoretical mod- 
els, like the macroscopic-microscopic (MM) calcula- 
tions to explain involve some prior knowledge of den- 



sities, single-particle potentials and other bulk proper- 
ties which may accumulate serious error in the largely 
extrapolated mass region of interest. Theypredict the 
magic shells at Z=114 and 1 84 fl H H H which 
could have surprisingly long hfe time even of the or- 
der of a million years [j, [M S S [l3l ■ Some other such 
predictions of shell-closure for the superheavy region 
within the relativistic and non-relativistic theories de- 



pend mostly on the force parameters 111, 121. 



Experimentally, till now, the quest for superheavy nu- 
clei has been dramatically rejuvenated in recent years 
owing to the emergence of hot and cold fusion reac- 
tions. In cold fusion reactions involving a doubly magic 
spherical targ et and a deformed projectile were used by 
GSI ill H [ll [13, [l3l to produce heavy elements 
upto Z = 110-112. In hot fusion evaporation reac- 
tions with a deformed transuranium target and a dou- 
bly magic spherical projectile were used in the syn- 
thesis of superheavy nuclei Z = 112-118 at Dubna 
[18, 19, 20, 21, 22, 23, 24J. At the production time of 
Z-\\2 nucleus at GSI the fusion cross section was ex- 
tremely small (Ipb), which led to the conclusion that 
reaching still heavier elements will be very difficult. At 
this time, the emergence of hot fusion reactions using 
*^Ca projectiles at Dubna has dramatically changed the 
situation and nuclei with Z = 114-118 were synthe- 
sized and also observed their a-decay as well as termi- 
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nating spontaneous fission events. It is observed that 
Z - 115-117 nuclei have long a-decay chains contrast 
to the short chains of Z = 114-118. Moreover, the 
life times of the superheavy nuclei with Z - 110-112 
are in milliseconds and microseconds whereas the life 
time of Z = 114-118 up to 30 s. This pronounced in- 
crease in life times for these heavier nuclei has provided 
great encouragement to search the magic number some- 
where beyond Z =1 14. Moreover, it is also an interest- 
ing and important question for the recent experimental 
discovery |25, 2^ 27 1 say chemical method of Z = 122 
from the natural 211,213,217,218'pjj which have long lived 
superdeformed (SD) and/ or hyperdeformed (HD) iso- 
meric states 16 to 22 orders of magnitude longer than 
their corresponding ground-state (half-life of ^"^^122 is 
fi/2 > 10^ years). 
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Figure 1: The proton average pairing gap Ap for Z= 11 2- 126 with 
N= 162-220 and Z= 112-130 with N= 162-260. 



In this letter, our aim is to look for the next double 
closed nucleus beyond ^o^p], vvhich may be a possible 
candidate for the experimentalists to look for For this, 
we have used two well-defined but distinct approaches 
(i) non-relativistic Skryme-Hartree-Fock (SHF) with 
FITZ, SIII, SkMP and SLy4 interactions S H (ii) 
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Figure 2: Same as FIG. 1 but for neutron average pairing gap A„. 



Relativistic Mean Field (RMF) formalism 1301 |31t] with 
NL3, Gl, G2 and NL-Z2 parameter sets. These mod- 
els have been successfully applied in the description 
of nuclear structure phenomena both in y6-stable and 
yS-unstable regions throughout the periodic chart. The 
constant strength scheme is adopted to take care of pair- 
ing correlation Il32i1 and evaluated the pairing gaps A„ 
and Ap for neutron and proton respectively from the 
celebrity BCS equations 113 311 . 

We scanned a wide range of nuclei starting from the 
proton-rich to the neutron-rich region in the superheavy 
valley (Z=l 12 to Z=130). It is well understood and set- 
tled that the properties of a magic number for a nuclear 
system has the following characteristics: 

• The average pairing gap for proton Ap and neutron 
A„ at the magic number is minimum. 

• The binding energy per particle is maximum com- 
pared to the neighboring one, i.e. there must be a 
sudden decrease (jump) in two neutron (or two 
proton) separation energy 52,, just after the magic 
number in an isotopic or isotonic chain. 
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Figure 3: The two neutron separation energy S2n and the shell correc- 
tion energy E shell for Z=l 12-126 and A'= 162-220 in the framework of 
SHF theory 



• At the magic number, the shell correction energy 
Esheii is maximum negative. In other words, a pro- 
nounced energy gap in the single-particle levels 
e„,p appears at the magic number 

We focus on the shell closure properties in the super- 
heavy valley based on the above three important observ- 
ables and identify the magic proton and neutron num- 
bers. 

The average pairing gap for proton Ap and for neutron 
A„ are the representative of strength of the pairing cor- 
relations. The curves for are displayed in FIG. 1 ob- 
tained by SHF and RMF with FITZ, SIII, SLy4, SkMP 
and NL3,NL-Z2, Gl, G2 force parameterizations. If we 
investigate the figure carefully, it is clear that the value 
of Ap almost zero for the whole Z=120 isotopic chain in 
both the theorical approaches. A similar A,, is observed 
for few cases of Z=124 and Z=l 14 isotopes. 

To predict the corresponding neutron shell closure of 
the magic Z=120, we have estimated the neutron pair- 
ing gap A„ for all elements Z=l 12-130 with their cor- 
responding isotopic chain. As a result of this, the cal- 



Table 1: Single-particle levels for neutron 6„ (MeV) for 120 in 
SHF(SLy4 and FITZ) and 3°" 120 in RMF (NL3 and G2). 



Orbit 


neutron (e^) 




SLy4 


FITZ 


NL3 


G2 


m 


-38.6 


-34.6 


-39.8 


-38.8 


pill 


-34.8 


-31.1 


-36.3 


-35.1 


1 11 

pl/2 


-34.6 


-31.0 


-36.1 


-34.8 




-29.9 


-26.6 


-31.4 


-30.2 




-29.2 


-26.1 


30.7 


-29.3 


s ' 


-26.2 


-23.1 


-26.3 


-26.1 


fill 


-24.2 


-21.3 


-25.7 


-24.5 


f5/2 


-22.7 


-20.2 


-24.2 


-22.8 


pill 


-19.1 


-16.5 


-19.8 


-19.0 


111 
p ' 


-18.9 


-16.3 


-19.7 


-18.7 


g9l2 


-17.9 


-15.3 


-19.3 


-18.1 


g7l2 


-15.3 


-13.4 


-17 


-15.4 




-11.9 


-9.5 


-12.9 


-11.6 




-11.1 


-8.8 


-12.5 


-11.3 




-10.9 


-8.7 


-12.3 


-10.7 




-9.8 


-7.2 


-10.2 


-9.3 




-7.3 


-6.0 


-9.3 


-7.5 


y7/2 


-4.5 


-2.4 


-5.8 


-4.1 


•13/2 


-4.0 


-2.0 


-5.5 


-4.1 


y5/2 


-2.6 


-0.9 


-4.7 


-2.4 


pill 


-1.4 


0.4 


-2.6 


-0.8 



culated A„ for the whole atomic nuclei in the isotopic 
chains are displayed in FIG. 2. We obtained an arc like 
structure with vanishing A„ at A^= 182, 208 and N-112, 
184, 258 respectively for SHF and RMF of the consid- 
ered parameter sets. Further, the neutron pairing gap is 
found to be minimum among the isotopic chains point- 
ing towards the magic nature of Z =120. Therefore, all 
of these force parameters are directing Z — 120 as the 
next magic number after Z =82. 

As mentioned earlier, the binding energy per par- 
ticle (BE/A) is maximum for double closed nucleus 
compared to the neighbouring one. For example, the 
BE/A with SHF (FITZ set) for 300,302,304j20 are 7.046, 
7.048 and 7.044 MeV corresponding to N=180, 182 and 
184 respectively. Similarly with SLy4 these values are 
6.950, 6.952 and 6.933 MeV. This is reflected in the sud- 
den jump of S2n from a higher value to a lower one at 
the magic number in an isotopic chain. This lowering 
in two neutron separation energy is an acid test for shell 
closure investigation. FIG. 3 shows the S 2,, as a function 
of neutron number for all the isotopic chain of the con- 
sidered elements for both SHF and RMF formalisms. In 
spite of the complexity about single-particle and collec- 
tive properties of the nuclear interaction some simple 
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Table 2: Same as Table 1. but for neutron e„ (MeV). 



orbit 


proton (ep) 




SLy4 


FITZ 


NL3 


G2 


^1/2 


-58.0 


-50.7 


-55.4 


-54.0 


P3/2 


-53.7 


-47.4 


-51.8 


-50.2 


Pl/2 


-53.4 


-47.2 


-51.6 


-50.0 


dsil 


-48.0 


-42.9 


-46.7 


-45.1 


di/i 


-47.2 


-42.3 


-46.0 


-44.3 


Sl/2 


-43.8 


-39.2 


-41.0 


-40.3 


fi/2 


-41.5 


-37.5 


-40.6 


-39.1 


/5/2 


-39.9 


-36.4 


-39.3 


-37.6 


Pyi 


-36.0 


-32.8 


-34.6 


-33.5 


Pl/2 


-35.8 


-32.5 


-34.5 


-33.1 


89/2 


-34.2 


-31.5 


-33.9 


-32.5 


8l/2 


-31.7 


-29.6 


-31.8 


-30.0 


dill 


-28.0 


-26.2 


-27.8 


-26.3 


di/i 


-26.8 


-25.2 


-27.2 


-25.4 


hn/2 


-26.5 


-25.0 


-26.9 


-25.3 


Sl/2 


-25.1 


-24.1 


-24.8 


-23.3 


hg/2 


-22.7 


-22.2 


-23.8 


-21.8 


fi/2 


-19.8 


-19.2 


-20.5 


-18.7 


in/2 


-18.5 


-18.1 


-19.6 


-18.1 


fi/2 


-17.7 


-17.5 


-19.4 


-16.9 


Pi/2 


-16.5 


-15.9 


-16.9 


-14.9 


Pl/2 


-16.2 


-15.7 


-16.7 


-14.4 


'11/2 


-13.3 


-14.1 


-15.6 


-13.1 


§9/2 


-11.7 


-11.9 


-13.2 


-11.0 


;i5/2 


-10.3 


-10.9 


-12.1 


-10.5 


§7/2 


-8.8 


-9.6 


-11.5 


-8.6 


d5/2 


-8.0 


-8.5 


-9.5 


-7.2 


d3/2 


-7.0 


-7.7 


-9.2 


-6.6 


Sl/2 


-3.6 


-5.7 


-8.2 


-6.0 


713/2 






-7.3 


-4.3 



phenomenological facts emerge from the bulk proper- 
ties of the low-lying states in the even-even atomic nu- 
clei. The 5'2« energy is sensitive to this collective/single- 
particle inter play and provides sufficient information 
about the nuclear structure effects. From FIG. 3, we no- 
tice such effect, i.e., jump in two neutron separation en- 
ergy at N=182 and 208 with SHF. However, such jumps 
are found at N=172, 184, 258 in RMF calculations con- 
firming the shell closure properties of the neutron num- 
bers. 

The shell correction energy Eshdi is a key quantity 
to determine the shell closure of nucleon. This con- 



cept was introduced by Strutinski 113411 in liquid-drop 
model to take care of the shell effects. As a result, the 
whole scenario of liquid properties converted to shell 
structure which could explain the magic shell even in 
the frame-work of liquid-drop model. The magnitude 
of total (proton plus neutron) E shell energy is dictated 
by the level density around the Fermi level. A positive 
Esheii reduces the binding energy and a negative shell 
correction energy increases the stability of the nucleus. 
As a representative case, we have depicted our SHF re- 
sult of E shell in FIG. 3. It is clear from the figure the 
extra stability of 302,328 220. We find similar results of 
large negative shell energy for RMF calculation at neu- 
tron number 172, 184, 258. Such calculations for few 
cases are reported in Ref. 13511 . 

As a further confirmative test, the single-particle en- 
ergy levels for neutrons and protons e„ ^ are analyzed. 
The calculated e„^p are presented in Table I for ^°^120 
SHF(SLy4 and FITZ) and for ^°^\2Q RMF(NL3 and 
G2) as representative cases. From the Table, one can 
estimate the energy gaps Ae„_p for neutron and proton 
orbits. For example, in ^"^120 (FITZ), the gap Ae„ - 
e„(M3/2) - e„(4.?i/2) at N=182 is 1.977 MeV which is 
comparable with 1.898 MeV of the known magic gap 
at N=50 for the same nucleus. For the proton case at 
Z=120, we get Aep = 6^(2/5/2) - fp(3p3/2) = 1-340 
MeV which is again a large gap of similar order (Ae^ - 
ep(l§9/2) - ep{^gi/2) = 1-862 MeV) forN=50. Almost 
identical behaviour is noticed with other SHF and RMF 
(at N=184) calculations, irrespective of parameter used, 
confirming Z=120 as a clear magic number. 

In summary, we have analyzed the pairing gap Ap and 
A„, two-neutron separation energy S2n, shell correction 
energy Esheii and single-particle energy for the whole 
Z=112-130 region covering the proton-rich to neutron- 
rich isotopes. To our knowledge, this is one of the first 
such extensive and rigorous calculation in both SHF 
and RMF models using a large number of parameter 
sets. The recently developed effective field theory mo- 
tivated relativistic mean field forces Gl and G2 are also 
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involved. Although the results depend slightly on the 
forces used, the general set of magic numbers beyond 
^"'^Pb are Z=120 and N=172, 182/184, 208 and 258. 
The highly discussed proton magic number Z - 1 14 in 
the past (last four decades) is found to be feebly magic 
in nature. 

We thank Profs. L. Satpathy, C.R. Praharaj and 
K. Kundu for discussions and a careful reading of the 
manuscript. One of the authors (MB) thank Institute of 
Physics for hospitality. 
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